Biaxial nematic and smectic phases of parallel particles with different cross sections 
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We have calculated the phase diagrams of one-component fluids made of five types of biaxial 
particles differing in their cross sections. The orientation of the principal particle axis is fixed in 
space, while the second axis is allowed to freely rotate. We have constructed a free-energy density 
functional based on fundamental-measure theory to study the relative stability of nematic and 
smectic phases with uniaxial, biaxial and tetratic symmetries. Minimization of the density functional 
allows us to study the phase behavior of the biaxial particles as a function of the cross-section 
geometry. For low values of the aspect ratio of the particle cross section, we obtain smectic phases 
with tetratic symmetry, although metastable with respect to the crystal, as our MC simulation 
study indicates. For large particle aspect ratios and in analogy with previous work [Phys. Chem. 
Chem. Phys. 5, 3700 (2003)], we have found a four-phase point where four spinodals, corresponding 
to phase transitions between phases with different symmetries, meet together. The location of this 
point is quite sensitive to particle cross section, which suggests that optimizing the particle geometry 
could be a useful criterion in the design of colloidal particles that can exhibit an increased stability 
of the biaxial nematic phase with respect to other competing phases with spatial order. 



I. INTRODUCTION 

The stability of thcrmotropic biaxial nematic phases 
0, 0] has been the subject of many studies in the last 
three decades, the driving force for this effort being their 
potential for use in display devices and other applications 
% Theoretical 0, H 1, Q and simulation @, 0, |, EH 
studies of a number of one- and two-component model 
systems have indicated that biaxiality is indeed possible 
in bulk nematic phases, but experimental identification 
is difficult and, in fact, very few positive reports have 
appeared in the literature [2|. 

Biaxial phases are a theoretical possibility in fluids 
made of particles that deviate from a cylindrical shape 
[l[, something that all molecules do to a larger or lesser 
degree. However, a biaxial particle geometry does not 
necessarily lead to a biaxial nematic phase: the degree of 
biaxiality has to be large enough, and all interactions 
have to be tuned optimally. Some early claims for a 
clear identification of a biaxial phase [ll|, EH, EH were not 
sufficiently substantiated [lij], and it has not been until 
recently that biaxial nematics in one-component fluids 
have been identified unambiguously [TEl . EH (although 
with a very low value of the biaxial order parameter) in 
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fluids with rigid bent-core (V-shaped) molecules. Stud- 
ies on V-shaped molecules, based on simple statistical 
models [l?], EH and computer simulations [IH , have pro- 
vided qualitative theoretical support for the experimental 
observations. However, it is difficult to isolate molecular 
shape as the crucial ingredient that causes biaxiality in 
low-mass molecular fluids, since specific interactions may 
play a role. 

The model systems on which we focus in the present 
paper are hard particles, which can be thought of as ide- 
alised (and quite faithful) representations of particle in- 
teractions in colloidal (rather than molecular) fluids con- 
sisting of colloidal particles with a shape anisotropy. In 
fact, novel methods to synthesise colloidal metallic parti- 
cles with a variety of shapes, from ellipsoids 20] and par- 
allelepipeds to rhombohedra and tetrapods [2ll 1221 . [23| , 
are now available, which opens up new avenues for theo- 
retical exploration. Particularly interesting are particles 
with rectangular shape (parallelepipeds), since the pres- 
ence of sharp edges, corners and flat sides may be the 
source of new types of ordering. Metallic nanorods with 
a rectangular cross section have been synthesised using 
different materials [U US HH, and various theoretical 
studies using different techniques and approximations in- 
deed point to complex phase behaviour [27], [H, [25^ [3(| 
for particles with square cross section. 

Biaxial phases in mixtures are also a possibility, as 
some theoretical [HI, HH, HH, H3] and simulation [1H stud- 
ies demonstrate; a positive identification on an experi- 
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mental rod-plate system also exists (36j . In mixtures, dif- 
ficulties are associated with competition between biaxial 
nematic ordering and nematic demixing transitions [37j ] - 
The stabilisation of phases with partial (smectic phase) 
or complete positional order (i.e. freezing) is certainly 
an effect that competes strongly with the formation of a 
stable biaxial nematic phase in one-component fluids, so 
that the window of particle biaxialities where a biaxial 
nematic phase can exist is predicted to be very narrow 
(1^ . Hard nanorod models seem to be ideal systems to 
study these problems, since one can focus just on particle 
geometry and the effect this has on phase behaviour. 

In a previous paper, Vanakaras et al. [38| have pre- 
sented theoretical and simulation results for a fluid of par- 
allel hard particles with rectangular cross section of ar- 
bitrary transverse aspect ratio. The results of Vanakaras 
et al. indicate that indeed a biaxial nematic phase can 
be stabilised at the expense of the smectic phase when 
the rectangular cross-sectional aspect ratio is sufficiently 
large. Mixtures of these particles have a considerably 
broadened region of stability for the biaxial nematic 
phase with respect to the pure-fluid case, which of course 
implies that mixing is a general mechanism to stabilise 
the biaxial nematic phase in these systems. 

In the present paper we revisit and extend the type 
of particles studied by Vanakaras et al. but consider 
only one-component fluids, again using the approxima- 
tion that particle interactions are completely hard. Since 
the Onsager-type theory used by Vanakaras et al. should 
(by construction) only provide a gross picture of the 
phase equilibria, we propose a sophisticated density- 
functional theory that overcomes some of the deficiencies 
of Onsager-type theories. Although our theory can be 
formulated for a general mixture, we particularize here 
to one-component fluids of particles whose main molec- 
ular axes are assumed to point along a specified direction 
(nematic director) but that possess a general cross sec- 
tion, characterised by a second molecular axis, that can 
freely rotate in the plane perpendicular to the nematic 
director. The theory is used, subject to some simplify- 
ing assumptions, to study the stabilisation of nematic 
and smectic mesophases for a number of particle geome- 
tries having different cross-sectional areas, such as the 
rectangular and elliptical, among others. Our theoretical 
scheme can therefore assess the relative stability between 
uniaxial and biaxial nematic, and uniaxial and biaxial 
smectic phases. 

Our proposal, based on the different phase diagrams 
presented, is that, by optimising the particle cross sec- 
tion (i.e. considering a wider range of geometries, not 
necessarily rectangular), one can also improve the stabil- 
ity range of the biaxial nematic phase in one-component 
fluids made of colloidal nanoparticles. A universal (i.e. 
independent of the type of particle) feature of the phase 
diagrams obtained is that the two (uniaxial and biax- 
ial) nematic and the corresponding smectic phases (four 
altogether) meet at a 'four-phase point', also observed 
by Vanakaras et al. for rectangular cross sections, re- 



sulting from the convergence of the corresponding four 
second-order transition lines separating pairs of phases. 
The location of this point, which is a stability bound- 
ary for the biaxial nematic phase, depends on the parti- 
cle geometry and this suggests a mechanism to enhance 
this stability. This result may be relevant for the de- 
sign and synthesis of colloidal particles exhibiting biaxial 
phases. An additional prediction of our theory is that, 
for particles with rectangular section and low transverse 
aspect ratio, a further, tetratic smectic phase, possess- 
ing four-fold symmetry in the transverse plane, appears 
in the phase diagram, albeit in metastable form. This 
phase is reminiscent of the corresponding tetratic ne- 
matic phase observed in two-dimensional fluids of hard 
rectangles [1^, H3, E[ • In fact, the topology of the phase 
diagram in the case of rectangular areas seems to have 
the same limit as the corresponding hard-rectangle fluid 
in two dimensions at high packing fractions, a result that 
reflects the dimensional crossover property exhibited by 
the density functional. 

The paper is arranged as follows. In the following sec- 
tion we introduce the density-functional theory, with rel- 
evant details on the numerical implementation relegated 
to the appendix. Sec. IIIII contains the results for all the 
particle geometries considered, separating the rectangu- 
lar geometry from the rest, but stressing the differences 
and similarities in phase behaviour. We end with some 
conclusions and perspectives for future work. 



II. FUNDAMENTAL-MEASURE FUNCTIONAL 
FOR PARALLEL PARTICLES 

The system to be studied consists of a fluid of hard 
biaxial particles with characteristic lengths L (parallel 
to the z axis), a\ and a 2. (both perpendicular to the z 
axis). The cross-section (transverse) area of the particle 
does not vary along the long axis of the particle, which 
is taken to lie along the z axis (primary nematic direc- 
tor); particles can otherwise freely rotate about this axis. 
Thus the fluid is described in terms of the density pro- 
file p{z,4>), with </> the azimuthal angle of the particle 
second axis with respect to a fixed direction in the trans- 
verse plane (secondary nematic director). Since the long 
axis is fixed, a trivial scaling along this direction follows, 
and the phase behaviour is not going to depend on L; 
however, we expect a strong dependence on the particle 
transverse aspect ratio k, = eri/o^ (cti being the larger 
size, along the particle second axis). 

We have chosen five different transverse sections, hav- 
ing symmetries as shown in Fig. [T] rectangle (R), semi- 
discorectangle (SDR, consisting of a rectangle capped 
with only one semi-disc in one of their ends), discorect- 
angle (DR, obtained from the previous one by adding 
another semi-disc at the other end), ellipse (E), and del- 
toid (D, composed of an isosceles triangle and its reflec- 
tion through its common base) . Note that both the latter 
and the rectangular geometries degenerate into a square 
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FIG. 1: (Color online). Sketch of particle transverse sections 
used in the phase-diagram calculations. The limiting case 
k = 1 is sketched, and for SDR the limiting case k = 1/2 is 
also shown. 



when k = 1, while ellipses and disco-rectangles both have 
a disc as a limit. Finally, the minimum value of aspect 
ratio for SDR is k = 1/2, corresponding to a semi-circle 
(see Fig. [J). 

In the following we obtain a fundamental measure den- 
sity functional for the smectic phase of a fluid composed 
of biaxial particles with a given cross section (in the sense 
explained above) in the approximation that the particle 
long axes are parallel. For this purpose we will adopt the 
projecting procedure to construct a functional for such a 
three-dimensional fluid starting from the corresponding 
two-dimensional functional for particles with the same 
section. Details of this formalism are given in Refs. [42j 
and H|. 

Since we consider the smectic as the less symmetric 
liquid-crystalline phase in the present calculations, the 
corresponding two-dimensional particles are not assumed 
to exhibit any spatial ordering in the xy (transverse) 
plane and, therefore, the excess part of the free-energy 
density can be constructed from scaled-particle theory 
(SPT) [the uniform limit of most fundamental-measure 
functionals] by (i) making the density profiles depend 
only on the z-coordinate, and (ii) defining two weighted 
densities, namely the local packing fraction, 



rz+L/2 

r)(z) = a p(z')dz', 

Jz-L/2 



(1) 



where p(z) = J dtfip(z, fa) and a is the cross-section area 
of the particle, and the two-particle weighted density 

pz+L/2 pz+L/2 p27T p2ir 

N 2 (z) = dzi dz 2 / d(f>i / d(j) 2 

Jz-L/2 Jz-L/2 JO JO 

x p{z 1 ,fa)p{z 2l (2) 

where fa 2 = 4>i—4>2 and we defined A(0) = A oxc (0)/2— a, 
with A cxc (<fi) the excluded area between two cross sec- 
tions. Thus, we obtain the following expression for the 



generating function of the excess free-energy density [4 

N 2 (z) 



$(«) = _2W] n [i_q( z )] 



T]{z)' 



(3) 



Using the projecting procedure [42j, |43| , the excess part 
of the three-dimensional free-energy density can be ob- 
tained using the formula 



(4) 



$(z) = A$ (z ), 

resulting in 

$(z) = -n(z) In [1 - r)(z)] + 

an(z)iV 2 (z) 
+ [I -viz)] 2 ' 

where the one-body weighted density n(z) is defined as 



n(z)n(z) + Nxjz) 
1 - T](z) 



(5) 



n(z) 



P z 



(6) 



while a new two-particle weighted density is obtained as 

l-Z + L/2 r27T r2TT 

Ni[z) = / dzi / d(j)i I d4>2p{zi,<j>i) 

Jz-L/2 Jo Jo 

x [p(z - L/2, fa) + p{z + L/2, fa)} A{fa 2 ). (7) 

The excess part of the free-energy functional for the 
smectic phase is then obtained as 



V 



dz$(z), 



(8) 



with Q(z) given by Eqn. ((5]), d being the smectic period. 
The ideal part is given by 



V 



dz 



dMz,0)[lnp(z,^)V-l], (9) 



with V the thermal volume. This ends the descrip- 
tion of the theoretical tools that will be used to study 
the phase behavior of biaxial particles with different 
geometries. The equilibrium state of the system fol- 
lows from minimisation of the total free energy density 
T = (J~id + ^cx) /V. The minimisation will be per- 
formed numerically, using a variational procedure, and 
adopting two different approximations for the parameter- 
izations of the smectic density distribution (both of which 
obviously contain the correct one-particle distribution of 
the higher-symmetry phases). These parameterizations 
are described in detail in Appendix [XI 



III. RESULTS 

This section is divided into two parts. In Sec. IIII Al 
we describe in detail the results obtained for a system of 
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hard biaxial parallelepipeds (rectangular cross section), 
while in Sec. IIIIBI we present the phase diagrams ob- 
tained for the other cross sections, stressing the most im- 
portant differences in phase behavior. The phases found 
in the phase diagrams and the notations used are: uniax- 
ial nematic (N), biaxial nematic (Nb), tetratic nematic 
(Nt), uniaxial smectic (Sm), biaxial smectic (Snas), and 
tetratic smectic (Suit)- As shown later, the main conclu- 
sion that can be drawn for this study is that the variation 
of the cross-sectional geometry has a dramatic impact on 
the relative stability of the Nb phase. This in turn sug- 
gests a relatively simple criterion, useful in the design of 
colloidal particles, to enhance the stability of the biaxial 
nematic phase with respect to non-uniform phases (such 
as the different smectic phases considered here). The 
underlying mechanism is alternative to that observed by 
Vanakaras et al. [Hj], where an increase of the Nb sta- 
bility follows by mixing two species with different sizes. 



A. Hard biaxial parallelepipeds 

Here the cross section is a rectangle with aspect ra- 
tio k — <J\/(J2- We begin with a comparison between 
our theoretical model and standard isobaric Monte Carlo 
(MC) simulations conducted on systems of N ~ 10 3 bi- 
axial parallelepipeds with their long axes parallel to the 
z axis. These systems require long times to equilibrate, 
so simulations in excess of 2 x 10 6 sweeps per particle are 
needed. The main goal here is to make the comparison 
at the level of the equation of state (EOS), and for this 
purpose we chose k = 2, i.e. a small value of the as- 
pect ratio (since this is a harder test for the theory than 
a large value). Minimisation of the functional was done 
using the decoupling approximation (see Sec. IA which 
gives exact results for the simple Sm symmetry, with no 
in-plane orientational ordering, but is only approximate 
when in-plane order builds up. As we cover the N and 
Sm phases, ranging from small to high densities (but al- 
ways below the Sm-SmT,B transition), we have used the 
Gaussian parameterization proposed in Sec. IA 21 

The resulting EOS, pertaining to the N and Sm 
branches, is shown in Fig. HJa), along with the simulation 
results. First, we comment on the latter. A compression 
run (filled circles) was performed from the low-density ne- 
matic phase. At r\ ~ 0.31 small-amplitude density waves 
began to develop, which gave rise to a fully developed 
smectic density distribution at r\ ~ 0.34 (the 'pretransi- 
tional' modulation is probably due to the small system 
size along the z direction) . The smectic structure exhibits 
no in-plane order of any kind (either translational or ori- 
entational), and thus can be identified as a standard Sm 
phase. However, at r\ > 0.37 some kind of translational 
order sets in, resulting in defected density distributions 
along z, still without any orientational order in the xy 
plane. These structures may result from a tendency of 
the system to develop crystalline order; the fairly low 
value of packing fraction at which this phenomenon oc- 



curs suggests that a plastic solid phase (with particles 
located at the nodes of a 3D lattice, but with their sec- 
ond axes randomly oriented) may be involved. An ex- 
pansion run (open circles in Fig. from a perfect 
(biaxial) crystal at high packing fraction also produces 
such defected structures and does not help to clarify the 
situation. However, the suggestion that a plastic solid 
might be stabilised is indirectly supported by theoretical 
calculations of the spinodal to a crystal phase (K), to be 
presented below. The main conclusion from the present 
simulation study is that the window of smectic stability, 
A77 ~ 0.03, is relatively small for moderate aspect ratios. 
Further study is required to obtain a more quantitative 
picture of the phase behaviour of this system in the crys- 
tal region. In any case, one can see that the comparison 
between theory and simulation is fair, as far as the value 
of the pressure is concerned. The location of the N-Sm 
spinodal point as predicted by simulations is rj ~ 0.34 
while the theory gives a value of r\ = 0.274; this is rea- 
sonably close to the simulation result. 

Our theory does not make any prediction on the tran- 
sition to a crystalline phase; in fact, the extension of the 
present model to include columnar or crystalline ordering 
is not a trivial task. Even if an approximate functional 
could be proposed, its numerical minimisation would re- 
quire huge numerical work. Therefore, in an effort to 
elucidate the system behaviour observed in the simula- 
tions for 77 > 0.37, we have performed a stability analysis 
in the framework of the restricted-orientation approxi- 
mation (Zwanzig approach) to estimate the location of 
the N-K phase transition. The approximation involves a 
constraint on the orientation of the particle second axis 
to lie parallel to either the x or y axes. In this context, 
a fundamental-measure density functional was obtained 
in Ref. [42j], which can be applied to study phases with 
any spatial symmetry, in particular the crystal. 

Fig. Hfb) is the phase diagram for hard parallelepipeds 
with small aspect ratios (between 1 and 3), as obtained 
from the Zwanzig approach. The continuous curves are 
the N-Sm and Sm-SmB spinodals, while the dashed curve 
is the spinodal instability in the N phase with respect to 
crystal fluctuations; these fluctuations are seen to corre- 
spond to a plastic solid, as suggested by the simulations. 
As can be seen from the figure, for the particular case k = 
2, the packing-fraction interval between the N-Sm and N- 
K spinodals is At? = r] K ~rj Sm = 0.3579-0.3013 = 0.0566, 
a result consistent with the simulations. However, this re- 
sult is to be taken with care, because the Sm-K transition 
is expected to be of first order (since both phases have 
different symmetries) and, consequently, the bifurcation 
analysis from the N to K phase is but a gross estimate 
for the location of this transition. All we can say for cer- 
tain is that the K phase bifurcates from the nematic at a 
packing fraction above (but close to) the N-Sm transition 
(although it is possibly metastablc within some density 
interval after bifurcation). These results also show that, 
for small aspect ratios, the Sme phase is unstable with 
respect to the K phase. 
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FIG. 2: (a) Theoretical equations of state of the N (dashed 
curve) and Sm (solid curve) phases of hard biaxial paral- 
lelepipeds with n = 2. Results from MC simulations are 
shown with symbols. Filled circles: compression run from 
nematic phase; open circles: expansion run from perfect bi- 
axial crystal at high pressure. Bifurcation points of the N-Sm 
second-order transition obtained by theory and simulation 
are indicated by filled and shaded arrows, respectively. The 
open arrow indicates the approximate Sm-crystal transition 
as obtained by simulation (see text), (b) Phase diagram ob- 
tained from the Zwanzig approximation for small values of 
k. Solid curves: N-Sm and Sm-Sms spinodals; dashed curve: 
N-K spinodal. 



For higher aspect ratios (for which we have not per- 
formed simulation studies), in particular for k — 15, we 
have found that the SniB phase is stabilized at low values 
of the packing fraction. To obtain the phase behavior for 
this particular case, we have taken advantage of the fact 
that the mean density is relatively small and, therefore, 
we have used the Fourier transform parametrization of 
Sec. IA 31 which gives a quasi-exact representation of the 
true density profile (numerical convergence is guaranteed 
in this regime of 7]). The results are shown in Fig. [3] (a), 
where the free-energy density of the N, Sm and Sme 
phases is plotted as a function of the packing fraction. 
The system exhibits a second-order N-Sm transition, a 
relative small window of Sm stability, and then a con- 
tinuous Sm-Sme transition, the latter being the stable 
phase for higher densities (up to the freezing transition) . 
Thus, the two-dimensional orientational ordering of the 
particle second axis appears in a continuous fashion with 
increasing density. In Fig. [3] (b) the smectic period is 
plotted as a function of packing fraction for both types 



FIG. 3: (a): Free energy density $* = fiTvjV vs. packing 
fraction r\ of the N (dotted curve), Sm (dashed curve) and 
Sme (solid curve) phases of particles with R cross section and 
with aspect ratio k = 15. (b): smectic period of Sm (dashed 
curve) and SniB (solid curve) phases. The symbols represent 
the N-Sm and Sm-Sme bifurcation points, respectively. 



of smectic phases. It is interesting to note that the pe- 
riod of the SniB phase decreases very slowly with density, 
compared with the corresponding behaviour in the Sm 
phase. 

With the aim to understand the non-uniform spatial 
and orientational correlations in the Sme phase, we have 
plotted in Fig. 2] the evolution of the density p(z) and 
the biaxial order parameter (see lA 3[) Ai(z) profiles with 
the mean packing fraction rj. While the density inho- 
mogeneities build up with packing fraction, the order 
parameter profile, although globally increasing, becomes 
flat as a function of z with increasing r\. Also, the profiles 
are out of phase, i.e. particles at smectic layers have a 
slightly lower orientational order than those situated at 
the interstitials. This is an interesting structural feature 
that points to a non-trivial coupling between layers via 
interstitial particles. However, this weak effect becomes 
less and less relevant as r\ increases, as revealed by the 
function Ai (z) becoming practically constant as a func- 
tion of z. The latter fact justifies a posteriori the use of 
the decoupling approximation to study the Sine phase 
and, probably, also other phases, such as the Smx phase 
that appears at higher densities and small values of k. 

To calculate the global phase diagram we have used two 
approximations: the decoupling approximation p{z, (/)) — 
p(z)h((f>) with a Gaussian parameterization of p(z) for 
k < 10, and the Fourier-transform parametrization of 
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FIG. 4: Evolution of the density profile (a) and biaxial order 
parameter (b) of the smectic phase of particles with R cross 
section and aspect ratio k = 15 for r\ = 0.16 + 0.02i (i — 
0, • ■ ■ ,4) and r] = 0.25 (from bottom to top). In (b) the case 
rj = 0.17 is also included. 



the complete function p(z, <j>) for k > 10. This choice 
is motivated by the dependence of the numerical crite- 
rion for convergence, in the minimization routines, on 
the value of r\. The phase diagram is plotted in Figs. [5] 
(a) and (b). For small aspect ratios, 1 < n < 2.618, we 
find a N-Sm transition at low densities, the Sm phase be- 
ing stable up to rj ~ 0.8, beyond which the fluid exhibits 
a continuous transition to a Smx phase. This transition 
was calculated using Eq. (|A11|) with i = 2. As described 
previously, the Smx phase consists of smectic layers in 
which the second particle axes (parallel to the smectic 
planes) point, with equal probability, along two mutu- 
ally perpendicular directions (secondary nematic direc- 
tors); within our approximation, the orientational distri- 
bution function fulfills the symmetry h{<f) = h((f> + 7r /2). 
This phase is sandwiched between the Sm and the Sme 
phases. Although we have not calculated numerically the 
location of the Smx-Sme transition, it can be approxi- 
mated, as noted in a previous work (40| , by the N-Sme 
spinodal extended to small values of k [the dashed curve 
of Fig. E (a)]. 

The Sm-Smr phase transition can be preempted by 
a transition to a crystalline phase with tetratic sym- 
metry. Nevertheless, we have shown the high density 
part of the phase diagram with only one-dimensional pe- 
riodic phases included. As we will see later, these re- 
sults are very useful for testing the performance of the 
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FIG. 5: Phase diagram for biaxial parallelepipeds (R cross 
section) calculated using the decoupling approximation with 
the Gaussian parametrization (a) and the Fourier expansion 
(b) for the density profile. The solid curves represent the 
continuous phase transitions between phases with different 
symmetries, as labelled in the figure. In (a) the dashed curves 
show the continuation of the Sm-Sme and Sm-Smi spinodals. 



present functional in the description of highly inhomo- 
geneous fluids. Further, for aspect ratios in the range 
2.618 < k < 18.101, the window of Sm stability (be- 
tween the N and the Sme phases) decreases with k, dis- 
appearing altogether at K = 18.101, the point where the 
N-Sm and Sm-SniB spinodals meet. This point was cal- 
culated using Eqns. (|A19p - (|A20p . with the inverse of the 
structure factor given by Eqn. (|A7|1 . For higher values 
of k the N phase exhibits a second-order transition to 
the Nb phase at a packing fraction calculated through 
Eqn. (j A19|) . On further increasing the density, there is a 
continuous transition between the Nb and SmB phases. 
The Nb-Shib and the N-Nb spinodals meet at the point 
mentioned in the introduction, which will be called four- 
phase point [38| since four different spinodals meet at the 
same point in the phase diagram. It is interesting to note 
that the transition between the Nb and SmB phases is 
reentrant in an interval of aspect ratios just below the 
four-phase point. This is a genuine prediction of our 
theory, since Onsager theory predicts a monotonic phase 
boundary between these two phases (3ST |. 

It should be noted that our prediction for the location 
of the four-phase point at k = 18.101, using our density- 
functional approximation, is to be contrasted with the 
value k ~ 15 reported in Ref. [38| . where a second-virial 
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Onsager theory was used instead. The MC simulations 
carried out in Ref. [38] seem to confirm this latter value. 
However, this apparent agreement should be taken with 
some caution. First, Onsager theory is known to give 
a poor picture of the N-Sm transition due to the mis- 
representation of density correlations. Also, excluded- 
volume effects underlying in-plane orientational ordering 
are probably not enough to give a quantitative descrip- 
tion, since higher-than-two-body correlations are known 
to be very important in two dimensions, and the problem 
at hand, once smectic layers have been formed, is quasi 
two-dimensional in nature. On the other hand, simula- 
tions of these systems are difficult, and large system-size 
effects are expected. The agreement found in Ref. [38| 
could be just fortuitous. Our approach, which includes 
higher-order correlations, should in principle give a more 
representative picture, but the situation is difficult to as- 
sess for lack of more extensive computer simulations and 
theoretical studies. For example, our theory gives only 
an approximate value for the third virial coefficient and, 
as shown in Ref. [44| , this coefficient is of the same order 
of magnitude as the second one for two-dimensional par- 
ticles (the particle cross sections) in the the limit k — > oo. 
Thus, for high aspect ratios our theory can deviate from 
the Onsager theory. A third-virial theory, including the 
exact virial coefficients up to the third order, is required 
to improve understanding of this issue. 
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FIG. 6: Phase diagram for particles with SDR cross section. 
Labels and lines are as in Fig. \S\ 



B. Other cross sections 

In this section we present the phase diagrams corre- 
sponding to the other particle cross sections described 
at the beginning of Sec. [Til As the symmetries of the 
different phases and the nature of their phase transitions 
were discussed in the preceding section, here we will con- 
centrate only on describing the differences between the 
phase diagrams of the different particle geometries. 

We begin by presenting the results corresponding to 
particles with SDR cross section (this is the only particle 
that does not possess head-tail symmetry). In Figs. [6] (a) 
and (b) we show the phase diagram. The location of the 
four-phase point, k = 17.f2f , is shifted to lower values 
of the aspect ratio, compared with the preceding case 
(rectangular cross section). Also, an important difference 
lies in the phase behaviour for small k and high packing 
fraction. To better visualise the difference, a zoom is 
shown in Figs. [3(a) and (b) in the region of Shit stability 
corresponding to R and SDR cross-sections, respectively. 

It is apparent from the figures that the stability of the 
SrriT phase (0.696 < k < 1.8f0) decreases by adding a 
semi-disc at one end of the rectangular section. This can 
be explained by the increasing excluded volume involved 
in the T-configuration (two particles in a perpendicular 
configuration) when this semidisc is added. It is also in- 
teresting to note that, for 0.5 < k < 0.696, a transition 
between the Sm and SniB phases again occurs, due to the 
fact that the rectangular part is so small compared with 



the total area of cross section that the T-configuration 
is not entropically favoured. Finally, in Figs. [7] (a) and 
(b) , the spinodals of the transitions between the isotropic 
phase (I) and the nematic and tetratic nematic phases 
of a strictly two-dimensional fluid of particles, with the 
same cross sections as analyzed here, are plotted. A rele- 
vant conclusion that can be drawn from the figure is the 
coalescence of these spinodals and those of the three- 
dimensional particles as packing fraction increases. This 
result confirms the dimensional cross-over property of 
the present density functional. For high i] the smectic 
phase can be considered as a collection of smectic layers 
where particles are perfectly located, and the effective 
interaction between particles located at different planes 
play a secondary role. 

The phase diagrams for particles with cross sections 
corresponding to DR and E geometries are shown in Figs. 
|5](a)-(b) and M (a)- (b), respectively. These two phase di- 
agrams have the common feature that the Smx phase 
is absent. The four-phase points are now located at 
k = 16.131 and n = 9.205, respectively. Finally, Figs. 
fTTJI (a)-(b) are the phase diagrams of particles with D 
section. Now at high packing fractions, and for aspect 
ratios in the range 1 < k < 1.430, the Sm phase exhibits 
a second-order transition to a Smx phase. However, the 
spinodal for this transition [^(k)] is an increasing function 
of k, which can be explained by the change of particle ge- 
ometry with k. The deltoid in the limit k = 1 coincides 
with a square; departure from this limit by increasing 
k involves a change in the angle between the adjacent 
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FIG. 7: Phase diagrams for R (a) and SDR (b) particles in the 
neighborhood of the Shit phase. All transition lines for three- 
dimensional fluids are shown with solid curves, while those 
corresponding to two-dimensional !4Qj with the same particle 
transverse sections are shown with dotted curves. For 2D sys- 
tems one needs to make the substitution Sm — > I (isotropic), 
Sme — > N (uniaxial nematic), and SmT — + Nt (tetratic ne- 
matic) . 



sides of the deltoid from its initial value of 90° and, as 
a consequence, the T-configuration of a pair of particles 
(which is characteristic of the tetratic symmetry) is less 
favoured. 

An additional feature that depends on the particle ge- 
ometry is the occurrence of reentrant behaviour in the 
transition between the Nb and Sme phases. This be- 
haviour is clearly associated with the rectangular nature 
of the particle shape, as it only appears in the phase 
diagrams of R and SDR particles, and much more pro- 
nounced in the former. Since both phases exhibit biaxial 
order, and the reentrant transition involves nematic to 
smectic ordering (or vice versa) , the effect is the result of 
a nontrivial coupling between spatial ordering along the 
z direction and angular ordering in the transverse plane. 

To explain the evolution of the four-phase points with 
the change in particle geometry we resort to Fig. 1111 
where the ratio between the coefficients — A* 2 and Aq is 
plotted for different cross sections. This ratio is a mea- 
sure of the relative reduction in excluded volume, or rel- 
ative gain in free volume, when the particles are orien- 
tationally ordered along the nematic director. As can 
be seen from the figure, this gain increases by modify- 
ing the particle sections in the sequence R, SDR, DR. E 




FIG. 8: Phase diagram of particles with DR cross section. 
Labels and lines are as in Fig. [3] 
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bottom to top). The location of the four-phase points are 
also shown (filled circles). The dotted curve represents the 
asymptote k — > oo for all the curves shown. 



and D. This in turn explains the sequence found in the 
location of the four-phase point, namely the Nb phase 
appears, for the first time, for the particle geometry that 
maximizes the gain in free volume associated to the ori- 
entational ordering of the second particle axis. 



IV. CONCLUSIONS 



In the present article we have analysed the phase be- 
haviour of models of particles that exhibit biaxial liquid 
crystalline order, with an emphasis on biaxial nematic 
phases. The models consist of hard particles with their 
principal axes parallel to each other, and such that the 
cross section along this common axis is constant while 
the secondary axis associated with this cross section can 
otherwise rotate freely in the plane. The positive identi- 
fication of a biaxial nematic phase requires that its sta- 
bility be compared with that of competing phases with 
spatial order, such as the smectic phase. In order to 
incorporate these phases into the theoretical scheme, a 
proper treatment of correlations has to be done. Since 
Onsager-type theories present severe defficiencies in this 
respect, we have developed a density functional, based 
on fundamental-measure theory, which makes a more 
appropriate treatment of correlations at high densities. 
The theory has been applied to study nematic and smec- 
tic phases with different orientational symmetries, such 
as the biaxial and tetratic symmetries, and the global 
phase diagrams, for particles with five different cross sec- 
tions. 

As our first and most important result, we have ob- 
tained the evolution in phase behavior with particle ge- 
ometry. For small aspect ratios and for R, SDR, DR 
and D-type sections, we have found a smectic phase with 
tetratic symmetry. The spinodals of the phase transitions 
between Sm and Sulb.t phases, at high packing fractions, 
are similar to those corresponding to phase transitions 
between isotropic and uniaxial or tetratic nematic phases 
of a strictly 2D fluid composed of particles with the same 
cross section. This result confirms the dimensional cross- 
over property of the functional. However, the Smx phase 
is preempted by the crystalline phase, as the MC simu- 
lations seem to show. All the phase diagrams for large 
values of k have a common feature, consisting in the pres- 
ence of a four-phase point at which the four spinodals 
corresponding to the second-order N-Sm, N-Nb, Sm- 
SniB and Nb-Sitib transitions meet. By studying the 
location of this point as a function of particle geometry, 
we have obtained a procedure to increase the Ns-phase 
stability, which might be useful in the design and syn- 
thesis of colloidal particles exhibiting a transition to this 
phase. In particular, the deltoid seems to be the cross 
section that favours the Nb stability most. Note that 
this idea is alternative to that proposed in Ref. [3&| . 
which consists in mixing two species with the same cross 
section but different particle lengths. 

The density functional proposed in this work can be 
used in a variety of situations, in particular, to study 
interfacial problems and the effect of confinement (for 
example in slit geometry) on the stability of different 
smectic phases of a fluid composed by biaxial particles. 
This study we leave for future work. 
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APPENDIX A: DENSITY PROFILE 
PARAMETERIZATIONS 

In IA II we describe the decoupling approximation and 
the corresponding expressions for free energy and struc- 
ture factor. In IA 21 a variational density profile, based 
on Gaussian trial functions, is proposed which, together 
with the decoupling approximation, allows to calculate 
the high density region of the phase diagrams. Finally, 
in IA3I a truncated Fourier expansion of the density pro- 
file is introduced; this parameterization is nearly exact 
for the description of smectic phases, but has the pit- 
fall that it can only be used for low values of the mean 
densities due to the poor numerical convergence of the 
minimization procedure. 



1. Decoupling approximation 

We adopt the usual decoupling approximation for the 
density profile: 



p{z,4>) = p(z)h(<j>), 



(Al) 



with h{4>) the orientational distribution function. As 
its name indicates, this approach decouples spatial and 
angular variables, which implies z-independent orienta- 
tional order in the smectic phase. In turn this means 
that the biaxial order parameters, defined by 



Ai = (cos(2i0)> = / d4>cos(2i<j>)h(<f>), (A2) 
Jo 

with i = l for uniaxial and i — 2 for two-dimensional 
tetratic symmetries, respectively, are constant within a 
smectic period. 

Inserting (|A1[) into ([5]), we obtain 



$(z) = -n(z) In [1 — r)(z)] + 
[1- V (z)} 2 ' 



n(z) V (z)(l + 2((A*))) 



1 - rj(z) 



(A3) 



where the double angular average 

{(A*))= dcf> 2 h(cj> 1 )h(h)A*((t> 1 -cp 2 ), 
Jo Jo 

(A4) 

has been defined. Also, the dimensionless quantity 
A*(<f>) = A(<f>)/a was introduced. Using the Fourier ex- 
pansion of the orientational distribution function, 



h(<f>) 



1 

2^ 



hk cos(k(j)) 



fe>i 



we obtain, for the double angular average, 



((A*)) = J2A* k h 



(A5) 



(A6) 



k>0 



with h = I. The Fourier coefficients A* k are given in 
Appendix [B] for all the particle sections studied (note 
that ^2n+i = f° r au tne geometries, except for SDR, 
the only one that breaks the head-tail symmetry; see 
discussion on the consequence of this in Sec. I A 3j) . 

The continuous N-Sm transition can be calculated 
from the divergence of the inverse structure factor 
S^ 1 (q,i]) = 1 — pc(q,r)), calculated from the Fourier 
transform of the direct correlation function c(q,r)); the 
latter is obtained from the second functional derivative 
of fiT with respect to the density profile. Within the 
decoupling approximation, this results in 

S- 1 ^) = 1 + 2yj 1 (q*) [2 + y + 2(1 + y) 2 ((A*))] 

+3i(q*/2)V [3 + 2y + 6(1 + y) 2 ((A*))] , (A7) 

with q* = qL, y = 77/(1 — 77), ji(x) = sin(x)/x, and where 
the double angular average ((A*)) is to be evaluated with 
h(<f>) = l/(27r), which gives the coefficient Aq. The equa- 
tion S^ 1 (q,rj) = 0, together with dS~ 1 (q,r])/dq = 0, 
must be solved for 77 and q at the absolute minimum of 
S-Hq,v)- 

2. Gaussian parametrization 



We adopt the following parameterized density profile: 

_l/2, 
v7T>' 



P 



(z) = pd (-) V2 ex P [- a ( z - kd f] > ( A8 ) 



i.e. a sum of normalized Gaussian peaks. This normal- 
ized form guarantees that <i -1 J Q dzp[z) = p, with p the 
mean smectic density. Insertion of (|A8[1 into ((T|) and (JfoJ) 

gives 



i(z) = ^7^{exp[-a(z + L/2-fcd) 2 ] 



exp [-a(z - L/2 - kdf] } 



r)(z) 



r/d 



^{erf [Va(z + L/2 - kd)} 



erf [y/a(z-L/2-kd)]}, 



(A9) 



(A10) 
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with erf(ir) the standard error function. For smectic sym- 
metry, and without orientational ordering parallel to the 
second nematic director, the expressions for ((A*)) with 
h((j)) — 1/2-7T (isotropic distribution function) are ana- 
lytic functions of the particle characteristic lengths and 
are provided in Appendix [B] for all the geometries con- 
sidered. We minimize the resulting free-energy density 
(3J-(a 1 d)/V with respect to the Gaussian parameter a 
and the smectic period d for a fixed mean packing fraction 
77 = paL. Varying 77 and repeating the above procedure, 
we obtain the free-energy branch for the Sm phase. The 
(continuous) nematic-smectic transition is located at the 
mean packing fraction value for which a ~ 0. Alterna- 
tively, this transition can be calculated from the diver- 
gence of the inverse structure factor, defined by (IA7I) . 

To calculate the second-order transitions between the 
Sm and the biaxial smectic (Snie) or tetratic smec- 
tic (Seit) phases, we have used a bifurcation analy- 
sis in which the orientational distribution function near 
the bifurcation point is approximated as h((f>) w [1 + 
hi cos(2i^>)]/(27r) (i = 1 and 2 for the uniaxial and 
tetratic symmetries, respectively). After insertion of this 
expression into the free-energy functional (3T obtained 
from the decoupling approximation (see Sec. IA lj) . we 
obtain the following expression for the free-energy dif- 
ference per particle (93 = [3T /N) between the Sm^ (or 
Shit) and Sm phases: 

h 2 h 2 ( AA* 

x /V(^H2-f)]\ (An) 

The spinodal curves (77 as a function of k) are then calcu- 
lated as the solution of the equation T(ri,a* ,d* ; k) = 
for rj, where a* and d* are those values obtained from the 
minimization of the free-energy density of the Sm phase 
f3T(a, d)/V with respect to the Gaussian parameter a 
and the smectic period d. 



3. Fourier parameterization and calculation of 
spinodals 

The density profile is now parameterized by a trun- 
cated Fourier expansion: 

K.M 

p(z,(f>) = — y~] s fcm cos(ofcz) cos(2to0), (A12) 

K,m>0 

where q — 2ixjd is the wave number. The latter, to- 
gether with the Fourier amplitudes s km , span the space of 
minimization variables. The zeroth-component Fourier 
amplitude is set equal to unity, i.e. s o = 1- Inserting 
the expression (|A12|) into the definitions of all the one- 



particle weighted densities, Eqns. ([I]) and ([6]), we obtain 
n(z) = p'^2s k ojo(qkL/2)cos(qkz), (A13) 

fe>0 

Vi z ) = V^2 s koji(qkL/2)cos(qkz), (A14) 

fc>0 

with jo(x) — cos a;. Also, we obtain the following ex- 
pressions for the two-particle weighted densities © and 
0: 

Ni(z) = 2pr] ^2 s k 1 nSk 2 nA* 2n j {qk 1 L/2)j 1 {qk 2 L/2) 

x cos(qfciz) cos(q£;2z), (A15) 
2 

N 2 (z) = — V s kin s k2n A 2n j 1 (qk 1 L/2)j 1 (qk 2 L/2) 
a L - J 

k\ ,/C2 ,71 

x cos(qkiz) cos(qk 2 z), (A16) 

where the coefficients A 2n for the different geometries 
are provided in Appendix [Bl The free energy is then 
minimised with respect to the Fourier amplitudes s km 
and to the wave vector q. In practice we needed about 
50 Fourier components with K = 10 and M = 5. 

In connection with the Fourier expansion for the den- 
sity (|A12[) and the corresponding expansions for the 
weighted densities (|A15p and (|A16|) . we must note that 
only coefficients of the excluded area with even index, 
A 2n , have been taken into account. This is justified for 
particles with head-tail symmetry, as A 2n+l = for 
these particles. However, for SDR particles, which do 
not exhibit this symmetry, one has A 2n+1 > 0. Since, 
for the nematic phase, or for the smectic phase in the 
framework of the decoupling approximation, we have 
{(A)) = ^2 k > A^h 2 , ^ ne free-energy minimization with 
respect to the Fourier amplitudes always gives h 2n +i = 0. 
However, if the coupling between spatial and orienta- 
tional degrees of freedom is properly taken into account, 
products of the form Sk im Sk2mA^ n , with odd m, do ap- 
pear in the expansions (|A15|) — (|A16[) for SDR particles. 
These coefficients could be negative, and in principle this 
could result in equilibrium values h 2n+ i ^ 0, i.e. in smec- 
tic phases with in-plane polar structure in their density 
profiles p(z, (f)). We expect terms of this type not to be 
dominant for large k, and therefore we have neglected 
these terms in the calculations for SDR particles, in the 
hope that the topology of the phase diagram near the 
four-phase point is not greatly affected by this approxi- 
mation. 

A measure of the local orientational order is given by 
the order parameters 

1 ^ 

A,-(z) = — - / d(j>p{z,4>)cos(2i(t>) 
P[z) Jo 

= TTT^y^SktCOsiqkz), (A17) 

2 pO) ^ 

with i = 1 for the uniaxial order and i — 2 for the tetratic 
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order; here the relation p(z) — p s^p cos(gfcz) should 

k 

be used. 

The second order N-Nb transition can be calculated 
using a simple bifurcation analysis of the free-energy dif- 
ference per particle between the phases, expressed as a 
truncated power series in the Fourier amplitudes hi (re- 
taining only the first term): 



as pointed out before, are also the equations, Eqns. 
(|A19I) - (|A20[) [obtained from the decoupling approxima- 
tion], whose simultaneous solution corresponds to the 
point where the N-Sm and N-Nb spinodal curves meet. 
This proves that this point is indeed a four-phase point. 



(A18) 



The non-trivial solution to the equation dA(p/dhi = 
gives 



y=yJl-(AA*) - 1, 
1 

n = l- 



(A19) 



The intersection between the spinodal of the N-Sm tran- 
sition, calculated using 



dq 



0, ((A*))=A* , (A20) 



[with the inverse structure factor S 1 (q, rj) given by (|A7[) ] 
and the spinodal of the N-Nb transition given ex- 

plicitly by Eq. (|A19[) ] can be found by substituting (|A19D 
into (|A7[) and solving (|A20[) for the variables q and k. Us- 
ing the Fourier parameterization approach the Sm-Sme 
transition is located at the value of 77 for which Ski ~ 0, 
Vfc, and finally the Nb- Shib bifurcation must fulfill the 
condition sxfe ~ Vfc. 

Let us now prove that the four-phase point must occur. 
In the neighborhood of this point, the leading order terms 
in the order-parameter expansion of the free-energy dif- 
ference between the Sulb and N phases from the Fourier 
parameterization approach [Eq. (|A12|) ] is 



= A(q, n; k)s 2 10 



-B(?7;«04l> 



(A21) 



where the coefficient A(q, 77; k) is proportional to the in- 
verse of the structure factor S~ 1 (q, rj) , given by Eq. (|A7|) , 



with ((A*)) = Aq, while B{t};k) is the coefficient of the 
quadratic term in the expansion of the free-energy dif- 
ference between the Nb and N phases with respect to the 
Fourier amplitudes hi, given by Eq. (|A18[) . Other terms 
in the expansion, depending on powers of sn and be- 
yond, are of higher order in magnitude. Thus, at this or- 
der of approximation, the smectic and orientational order 
parameters are decoupled. Minimization of (|A21[) with 
respect to S10 and sqi gives the set of coupled equations 



A(q*,r,;K) = 0, B(»j;k)=0, 



(A22) 



with q* the value at the absolute minimum of S~ 1 (q,rj) 
with respect to q. These equations should be solved to- 
gether for 77 and k to find the (unique) point in the rj— 
k plane where the N phase becomes unstable with re- 
spect to Siiib fluctuations. But note that Eqn. (|A22|) . 



APPENDIX B: EXPRESSIONS FOR THE 
COEFFICIENTS A* 2n 



In this section we provide the analytic expressions 
for the coefficients A\ n corresponding to all the particle 
transverse sections studied. 



1 



R : A* 2n = -{(k + k' 1 +2) S n 



l[ K+K - 1 + 2(-ir] { 4-^4\. ,Bi, 



SDR 



1 tt\ 1 (I 



2 8 



4n a - 1 , 

7T 2K 2 



k -\ 

2 V 8 



SnO 



1 / K+ ( - i) ;- i Vii^i, ( B2 ) 



2tt 



DR: A* = 



7T ( K — 1 



An A - 1 



( K -l) 2 (l-<5„ ) 



(An 2 



D: A* 2n 



( ~k 

2tt I K - 1 + - 

V 4 

(k + k- 1 ) 



(B3) 



2^no — cos(2n7) 



2(1 - <W 



An 2 - 1 



(B4) 



with 7 = arctanK 1 . Finally, the expression for the E 
geometry can be calculated only numerically from 



E : At 



(1 + S n0 ) f 71 
2tt 



d(f>cos(2n(f>) A* (</>), 



A*(cf>) = [BlW+aaW] £?[r(0)], 



(B5) 



with E(x) = J^ 2 dty/l — x sin 2 t the complete elliptic 
integral of the second kind and where there were defined 



Si((/>) = (k 2 - 1)| sine 



S-2 



(0) = v/^ + W, 



t(4>) 



\si{4>) + s 2 (4>)f 



(B6) 



For n = we have A* = Ak[E(\ - «r 2 )] 2 /7r 2 
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